Interplay between coarsening and nucleation in an Ising model with dipolar 

interactions 
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We study the dynamical behavior of a square lattice Ising model with exchange and dipolar 
interactions by means of Monte Carlo simulations. After a sudden quench to low temperatures we 
find that the system may undergo a coarsening process where stripe phases with different orientations 
compete or alternatively it can relax initially to a metastable nematic phase and then decay to the 
equilibrium stripe phase through nucleation. We measure the distribution of equilibration times for 
both processes and compute their relative probability of occurrence as a function of temperature and 
system size. This peculiar relaxation mechanism is due to the strong metastability of the nematic 
phase, which goes deep in the low temperature stripe phase. We also measure quasi-equilibrium 
autocorrelations in a wide range of temperatures. They show a distinct decay to a plateau that we 
identify as due to a finite fraction of frozen spins in the nematic phase. We find indications that 
the plateau is a finite size effect. Relaxation times as a function of temperature in the metastable 
region show super- Arrhenius behavior, suggesting a possible glassy behavior of the system at low 
temperatures. 
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I. INTRODUCTION 

It is well known that the usual long range order charac- 
teristic of a ferromagnetic phase dominated by exchange 
interactions can be destroyed by the presence of dipo- 
lar interactions inducing the formation of magnetic do- 
mainsi. In some cases, typically when perpendicular 
anisotropy is present, the dipolar interactions are anti- 
ferromagnetic in character and compete with the ferro- 
magnetic exchange interactions — >^ giving rise to antifer- 
romagnetic or stripe phases. These phases are common 
in ferromagnetic thin films with perpendicular anisotropy 
and have been the subject of intense experimental^i^i^i^ii^, 
theoretica l^i^°i^^ and numerica l^^i^^i^"^'^^'^^'^''^ work in the 
last 20 years. In the strong anisotropy Hmit an ultrathin 
film can be modeled by a system of magnetic dipoles 
on a square lattice in which the magnetic moments are 
oriented perpendicular to the plane of the lattice, with 
both nearest-neighbor ferromagnetic exchange interac- 
tions and long-range dipole-dipole interactions between 
moments. The thermodynamics of this system is ruled 
by the dimensionless Ising Hamiltonian: 

where S stands for the ratio between the exchange Jo > 
and the dipolar Jd > interactions parameters, i.e.. 



5 = Jo/Jd- The first sum runs over all pairs of near- 
est neighbor spins and the second one over all distinct 
pairs of spins of the lattice; is the distance, measured 
in crystal units, between sites i and j. The energy is mea- 
sured in units of Jd- The equihbrium phase diagram of 
this system has been extensibly studie d^^'^^'^^'^^ , while 
several dynamical properties at low temperatures were 
studied ii>i^ii2ii2£i2i. The threshold for the appearance 
of the stripe phase in this model is Sc — 0.425 ^^'^^ . For 
S > Sc the system presents a sequence of striped ground 
states, characterized by a constant width h, whose value 
increases exponentially with 6^^. 

Ini^ it was shown that in the range 1 < S < 3 the 
system presents a first order phase transition between 
a low temperature stripe phase and a high temperature 
tetragonal phase with broken translational and rotational 
symmetry. In a subsequent workii it was shown that 
for a narrow window around S — 2 the model shows 
an intermediate nematic phase, between the stripe and 
tetragonal ones, where the system has short range po- 
sitional order but long range orientational order. The 
nematic phase exists between two critical temperatures 
Ti < T < T2. For finite sizes these are actually pseudo 
critical temperatures Ti{L) and 12 (L), with finite en- 
ergy barriers between phases. The extrapolation to the 
thermodynamic Hmit gives ri(oo) ~ 0.772 and T2(oo) « 
0.797ii. For L - 50 (L = 48 - 56) we can estimate 
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Ti{L) and 72 (L) as the average between the positions of 
the minimum of the cumulant and the maximum of the 
specific heat, so Ti(50) « 0.78 and T2(50) « 0.81. For 
the range of system sizes considered {L ~ 50) the energy 
barriers between the nematic and the tetragonal phases 
are rather small, while the barriers between the stripe 
and the nematic phases are very largely. Hence, under a 
cooling from high temperature no metastable tetragonal 
states are observed, but we see a strong metastability in 
the nematic phase. This can be observed in Fig. 19 of 
Ref.— . Although metastability strongly resembles a first 
order phase transition, we were not able to determine the 
nature of the stripe-nematic transition unambiguously: 
while the nematic order parameter presents a jump at Ti 
and the free energy barriers seem to diverge in the ther- 
modynamic limit, the internal energy is continuous at the 
transition and the peak in the specific heat tends to satu- 
rated^. Further MC simulations showed that the nematic 
phase is present also for other values of (5 7^ 2 in differ- 
ent parts of the phase diagram, mainly around the phase 
borders between striped states of different widthslS-. 

Here we study the out of equilibrium dynamical prop- 
erties of the two dimensional Ising model with exchange 
and dipolar interactions described by the Hamiltonian ^ 
with 6 — 2. We performed quenches of the system from 
initial conditions at high temperatures directly to the 
stripe phase (T/ < Ti). Our main result is the identifica- 
tion of two kinds of processes: a slow coarsening towards 
a metastable nematic phase followed by nucleation of the 
stripe phase in a background of the metastable nematic 
and also direct coarsening of the stripe phase without an 
intermediate nucleation process (see Fig[T]). In each real- 
ization the system chooses one of both paths to the final 
stripe phase with probabilities that depend on the final 
temperature of the quench. We computed these prob- 
abilities from Monte Carlo simulations and then using 
arguments from homogeneous nucleation theory, finite 
size scaling analysis and known growth laws for critical 
dynamics we were able to characterize completely both 
kinds of processes in the whole temperature range below 
the stripe-nematic transition. 

We also made slow cooling experiments, similar to 
those already discussed inii and which show strong hys- 
teretic behavior in the energy. We confirm that hys- 
teresis is associated with a strong metastability of the 
nematic phase, while the tetragonal phase shows neg- 
ligible metastability in the energy under cooling. Fi- 
nally, (quasi-)equilibrium autocorrelation functions ob- 
tained for a whole range of temperatures between the 
tetragonal and the stripes phases show slow relaxation 
characterized by stretched exponential decay. In the 
stripe phase the relaxation times obtained from the au- 
tocorrelations exhibit a strong growth which can be well 
fitted by a Vogel-Fulcher-Tamman relation, similar to the 
behavior of fragile glass formers, with a relatively high ( 
To ~ 0.45) divergence temperature. 

In Section II we describe the quench experiments and 
the associated phenomenology. In Section III the slow 



cooling experiments and the behavior of the equilibrium 
correlations are discussed. In Section IV we present a 
discussion of results and brief conclusions. 



II. RELAXATION AFTER A QUENCH FROM 
HIGH TEMPERATURE: INTERPLAY BETWEEN 
COARSENING AND NUCLEATION 

We performed quenches from infinite temperature 
(completely random initial configuration) down to tem- 
peratures T < Ti. We also compared with the case 
where the initial configuration is taken from an equili- 
brated state at temperatures T ~ 1. The results were 
the same. 

The typical relaxation behavior can be appreciated in 
Figdl where we plotted the time evolution of the instan- 
taneous energy per spin E{t)/N for two different realiza- 
tions of the stochastic noise, together with typical spin 
configurations along the evolution. We see that the sys- 
tem can relax through two different types of mechanisms. 
In Fig[T^ the system relaxes first into a metastable ne- 
matic state, where it stays during a long time r, after 
which suddenly relaxes to the equilibrium striped state. 
In other words, there is a two-step relaxation: first a 
coarsening process, where domains of the nematic phase 
compete (horizontal and vertical stripe orientations) , fol- 
lowed by a nucleation of the stable phase (we will verify 
this quantitatively later). In Fig[TjD the system form do- 
mains of both the stable (striped) and metastable (ne- 
matic) phases that compete and relax directly to the 
equilibrium state through a coarsening process. Both 
types of relaxation processes can happen with certain 
probability which depends on the temperature. To ana- 
lyze quantitatively these processes we computed the equi- 
libration time probability distribution. 

The equilibration time r is defined as the time (in 
Monte Carlo Steps, MCS) that the system takes to reach 
the stable, crystal-like striped state after a quench. The 
criterium to determine r was the following. For each 
quench temperature T < Ti we calculated the equilib- 
rium energy per spin histogram along a large single MC 
path starting from the ground state. This histogram 
presents two gaussian-shaped peaks, each one centered 
at the averages Ustr{T) and Udis{T)r^. Between the 
two peaks the histogram presents a finite minimum (at 
least for a finite system), associated with the free en- 
ergy barrier that separates both phases. We also mea- 
sured the standard deviation of the low energy peak 
associated with the striped phase astr{T) and defined 
emin{T) = Ugtr + <^ str , verifying that it is well below 
the central minimum of the histogram; Cmm is almost 
independent of the system size for L > 32, so it can 
be easily calculated using small system sizes. Then for 
each quench we calculated the average energy per spin 
e{t) = {H) /N. When e{t) fell bellow e^in we stopped 
the simulation defining e(r) = Cmin (see Fig[l|). Repeat- 
ing this procedure we obtained the probability distribu- 
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FIG. 1: Time evolution of the energy per spin in a single MC run {S = 2, L=52 and T = 0.65) for two different realizations of 
the stochastic noise. Some snapshots of the spin configurations at different times are shown below the figures, (a) Competition 
between nematic domains (coarsening), that leads to a metastable nematic state, followed by an activated relaxation to the 
striped state (equilibration through nucleation); (b) Competition between nematic and striped domains (coarsening) that leads 
directly to equilibration. 



tion (normalized frequency) P{t) of the stochastic vari- 
able T for different values of T and L. 

The above mentioned relaxation mechanisms are re- 
flected in a characteristic two-peak structure in P(t), 
each one centered at typical values corresponding to well 
different time scales, as shown in Figsl2]and[3l This prob- 
ability distribution can be very well fltted using a super- 
position of two log-normal functions (see Figsl2] and [3]) 
P(lnT) = Pi(lnr) + P2(lnT) with 



p(lnT) = At cxp 



2 a? 



^ = l,2 (2) 



where Ai are normahzation constants. Using these fit- 
tings we can see in Fig. |4] that the characteristic time 
scales of both peaks increases as T decreases. 

These fittings also allowed us to estimate the averages 
(r)j = exp (^fii + / 2) , i = 1,2, associated with each 
process. Both averages increase with L, so we can use 
finite size scaling analysis to determine the nature of the 
associated processes. 
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FIG. 2: (Color online) Normalized frequency P(t) vs Inr 
for 5 — 2, T = 0.65 and different system sizes. The lines 
correspond to a fitting with a superposition of two log- normal 
distributions. Sample size for the histograms were ~ 2 x 10*. 
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FIG. 3: (Color online) Normalized frequency P(t) vs Inr 
for S = 2, T = 0.77 and different system sizes. The lines 
correspond to a fitting with a superposition of two log- normal 
distributions Eqs.(l2|). Sample sizes for the histograms run 
between ~ 4 x 10=* (L = 44) to ~ 5 x 10* {L = 16).) 




FIG. 4: Fitted probability distribution P(r) as a function of 
Inr and T for (5 = 2 and L = 24. 



A. Finite size scaling of the equilibration times 

Suppose that the relaxation is completely governed by 
a competition between ground state domains, that is, by 
a curvature driven relaxation of the domain walls. The 
excess of energy respect to the ground state is given by 
Se = jn, 7 being the "surface" tension per unit length 
of the domain walls and n the fraction of spins belonging 
to them; n is approximately given by n ~ bNd{t)l{t)/N, 
where Ndit) is the average number of domains, l{t) is the 
average linear size of the domains and & is a geometrical 
factor. Using Nd oc N/l'^{t) and grouping all these factors 
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FIG. 5: (Color online) (r)j^ vs. L from two log-normal fittings 
Eqs.© 



we get 5e ~ l^'^{t). For 5 — 2 the energy distribution of 
the striped phase is very narrow and lies very near the 
ground state, so the criterium e(T) — Cmin implies Se <C 
1, condition that is attained when 1{t) ^ L. Since for a 
heat bath dynamics with non-conserved order parameter 
(class A system) we can assume l{t) ~ t^/^ we expect the 
scaling (t) ~ L^. We verified that (r)^ ^ for a wide 
range of values of T and L (see an example in FigH]), 
thus confirming that Pi{t) is associated with a simple 
coarsening process. 

Now suppose that relaxation is governed by an homo- 
geneous nucleation process. According to classical nucle- 
ation theory (r) follows approximately an Arrhenius law 
(r) = To exp(A(L,r)/T), A(i,T) being the height of 
the free energy barrier to nucleation. In a finite two- 
dimensional system^^ A(Tm) ~ L, Tm being the melting 
temperature; in our case Tm = Ti. In an infinite sys- 
tem below the melting point, classical nucleation theory 
predicts that A(r) = A(oo, T) is given by the excess of 
free energy of a critical droplet with average linear size 
^(T) of the ordered (striped) phase inside the disordered 
(nematic) metastable state. Thus, if L £,{T) the free 
energy barrier is expected to be almost independent of 
the system size when T < Tm- But, if L < ^(T), the 
droplet never reaches the critical size and the excess of 
free energy is dominated by surface tension. Hence, the 
transition to the crystal phase will occur when the crystal 
droplet reaches the system size and therefore we expect 
A(L,r) - L and In (r) - L/T. Fig [6] shows that {t)^ 
follows such scaHng, suggesting its identification with the 
average nucleation time Tnuci- However, due to the pro- 
cedure we used to obtain it, {t)^ actually corresponds to 
the crystallization time, which is not necessarily equal to 
Tnuci^ because after the creation of a critical droplet it 
follows a growth process until the whole system reaches 
the crystal phase (for the rather small system sizes we 
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FIG. 6: (Color online) Estimated free energy barrier 
A{L,T) ?a Tin {{t)^) vs. L for different temperatures. The 
dashed lines correspond to a sigmoidal fitting. The inset 
shows the extrapolated free energy barrier A(r) = A(oo,T) 
vs. T. The error bars are shown only when larger than the 
symbol size. 



considered we can assume that there is only one critical 
droplet). Hence, {t)^ = Tnud + Tgrowth- However, we 
can expect Tgrowth ~ i", so Tnuci > Tgrowth, at least 
as far as L is large enough that we are in the scaling 
regime, but small enough so that the free energy bar- 
rier still depends on L. Under these conditions {t)^ 
give us a good approximation of Tnuci and we can esti- 
mate A(L,T) « T In (t)2 up to logarithmic corrections. 
Then, by fitting the curves of FigE] using a sigmoidal 
function A(L,T) = A(T)/ (1 + exp (i*(r) - L)/c(r)) 
(L*(T), A(T) and c(T) are fitting parameters) we can 
extrapolate A(T) = lim_L-^oo A(L, T). The extrapolated 
curve A(r) is shown in the inset of FigUl We see that 
A(T) grows strongly as T ^ , while it appears to 
present a minimum around T = 0.5, or at least to sat- 
urate in a constant value as T ^ 0. This is compatible 
with the results in^^ obtained from a different analysis 
where it was found a divergent free energy barrier at the 
transition temperature Ti in the thermodynamic limit. 
For quench temperatures T < 0.5 very large relaxation 
times make it very difficult to obtain reliable numerical 
results. For future analysis we list in Table U the approx- 
imated fitted values of 
T. 

We also calculated the relative probability of occur- 
rence of nucleation and coarsening processes 



t)^ for different values of L and 



nucl 



fZ,P2{\nT)dilnT) a^A^ 



Pn. 



J^^Pi{\nT)d{lnT) aiAi 



(3) 



in the range 0.5 < T < Ti. In Figl?] we see that (f>{T) 
decreases monotonously with T. It seems to diverge as 
T Tj" and L ^ oo, while becoming almost indepen- 
dent of L for T < 0.6, saturating in a constant value 
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FIG. 7: (Color online) Relative probability Eq.© vs. T for 
(5 = 2 and different values of L. 



4> ~ 0.5. This shows that nucleation dominates the relax- 
ation near the melting point, while coarsening becomes 
dominant for deep quenches, although there is always a 
finite and relatively high (~ 1/3) probability of nucle- 
ation. 



III. RELAXATION AFTER A SLOW COOLING 
A. Energy 

We first analyze the behavior of the average energy 
per spin during a slow cooling; that is, we first equili- 
brate the system at some temperature Tq > T2 and then 
decrease the temperature down to a final value T/ <C Ti 
according to the linear protocol T{t) =To — rt, where r 
is the cooling rate, t is in MCS and the average is taken 
over several runs between Tq and Tf. In Figl8]we show 
an example. We see that below some value of r the en- 
ergy curve becomes almost independent of r. In FiglH 
we compare the cooling curves with the corresponding 
heating curve; this last curve is calculated by perform- 
ing a Hnear heating from the ground state, starting at a 
very low temperature. We also compare the curves with 
the equilibrium values obtained in Ref?l^. These results 
verify that (for this range of system sizes) metastability 
in the tetragonal-nematic transition is negligible, while 
there is a strong metastability associated with the stripe- 
nematic transition. Notice that the metastable cooling 
curve has an infiection point around T = 0.5 and falls 
down; although it does not reach the equihbrium value 
for T < 0.5, this infiection point suggests the existence 
of a spinodal temperature around T = 0.5. Moreover, 
from Tabled] we see that for L 50 the nucleation time 
is Tnuci ~ 4 X 10^ MCS; then, for a cooling rate r = 10~^, 
the time elapsed from the instant at which the temper- 
ature r = Ti to the instant corresponding to a final 
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r = 0.60 


r = 0.65 


T = 0.70 


T = 0.75 


T = 0.77 


L = 32 


~ 1.3 X 10'' 


~ 5.4 X 10^ 


~ 2.8 X lO'^ 


~ 1.6 X lO'^ 


~ 2.0 X 10® 


L = 36 


~ 2.4 X 10*^ 


~ 1.1 X 10® 


~ 5.7 X 10^ 


~ 3.3 X 10® 


~ 4.7 X 10® 


L = 40 


~ 3.4 X 10*^ 


~ 2.0 X lO'* 


~ 1.1 X 10® 


~ 8.0 X 10® 


~ 1.0 X 10® 


L = 44 


~ 5.3 X lO'' 


~ 3.0 X 10® 


~ 2.0 X 10® 


~ 1.4 X 10® 


~ 1.6 X 10® 


L = 48 


~ 6.6 X 10*^ 


~ 4.4 X 10® 


~ 3.3 X 10® 


~ 2.1 X 10® 





TABLE I: (r)2 for different values of L and T. 




temperature around T — 0.5 will be t = AT/r ~ T„„ci- 
Therefore, for that range of temperatures we would ex- 
pect that in a large fraction of the realizations the system 
would have already equilibrated. Indeed, the results of 
FigdO] verify this assumption. Moreover, the tracking of 
typical spin configurations for those realizations were the 
system did not equilibrate below T = 0.5 shows that ac- 
tually the system is no longer in a nematic state; instead 
of that, it got stuck in another metastable state, with a 
mixture of stripes with widths h — 2 and h = 3. Actually 
the /i = 3 striped phase becomes metastable precisely for 
that range of temperatures^^. This behavior reinforces 
the hypothesis of the existence of a spinodal tempera- 
ture Tsp « 0.5. Also the behavior of the orientational 
order parameter observed in Fig. 19 of Refpi^ supports it. 
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B. Correlations 



FIG. 10: (Color online) Energy per spin cooling curves: indi- 
vidual realizations (5 = 2 and L = 48). 



We now analyze the two-times correlation function 



after a cooling from high temperature to a fixed tempera- 
ture T < T2. In this case we performed the cooling with 



a ladder protocol, that is, we reduced the temperature 
by steps of AT = 0.01, letting the system to equilibrate 
during 10'* MCS at each intermediate temperature. On 
the average, this corresponds to a cooling rate r — 10~^. 
Once we arrived to the measuring temperature (and af- 
ter another 10^ MCS) we set t = and we saved the 
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FIG. 11: (Color online) Two-times correlation function 
C{tn,,tw +t) as a function of t for 5 = 2, L = 48 and differ- 
ent waiting times t^, after a quasi-static cooling at T = 0.78. 
Sample size for average was 150 
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FIG. 12: (Color online) Two-times correlation function 
C{tn,, tw -I- t) as a function of t for 5 = 2, L = 48 and differ- 
ent waiting times tw, after a quasi-static cooling at T = 0.7. 
Sample size for average was 350 



initial configuration. The autocorrelations ([4]) were then 
measured as a function of t and and the whole curve 
was averaged over different realizations of the stochastic 
noise, starting always from the same initial configuration 
(this is to save CPU time). Some check repeating the 
calculation starting from different initial configurations 
(all obtained with the same cooHng protocol) showed no 
difference with the previous one. 

First of all we analyzed a cooling down to a tempera- 
ture Ti(L) < T < 12 (L), which corresponds to the stable 
region of the nematic phase. From Fig[lT] we see that 
C{tuj,tw + t) = C{t), as expected in a stable phase. We 
also noticed that the correlations do not decay to zero 
but to a plateau Cpiateau ~ 0.37. We will discuss the 
meaning of that plateau later. 

In Fig[T2]we show the same calculation, but for a tem- 
perature T < Ti{L). Both t and were chosen such t 
t + tw<^ Tnuci- We see that again C{tw,tw +t)=C 
consistently with the expected quasi-stationary natur 
a metastable state. We also see that the correlation ; 
displays the characteristic plateau observed in the st£ 
nematic state. 

In FiglH we show C{t) for T = 0.7 and L = 48 i 
very long run with t up to 2 x 10'' MCS. Comparing v 
table [J we see that the second decay in the correlal 
coincides with the average nucleation time. 

In FigdHwe show C{t) for L = 36 and different coo' 
temperatures. We see that the plateau appears for 
temperature, but becomes higher in the nematic cast 

What is the origin of the plateau? In FigHH we si 
Cpiateau as a function of L for different temperatures, 
temperatures in the tetragonal liquid region, we see t 
for large enough sizes the plateau decays as 1/y/N. 
expected in a disordered state. However, for small Siz.co 
the plateau is independent of L and the crossover size 
increases as the temperature approaches T2. Once we 
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FIG. 13: Correlation function C{t) as a function of ttor 5 — 2, 
Z/ = 48 and different waiting times t^, after a quasi-static 
cooling at T = 0.7. Sample size for average was 80 
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FIG. 14: (Color online) Correlation function C(t) as a func- 
tion of t for 5 — 2 and L — 36 after a quasi-static cooling at 
different final temperatures. 
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FIG. 15: (Color online) Correlation function plateau Cpiat 
as a function of L for 5 = 2 and different temperatures. 



are in the nematic phase (either stable or metastable) 
the plateau for small sizes becomes independent of T 
(at least for a range of temperatures where the relative 
probability of nucleation is high, see Fig. [7l and it is 
roughly equal to 1/3. Since in a stationary state the un- 
connected correlations we are considering here behave as 
C{t) ^ X^i i^i) 1 this result shows that in the nematic 
state one third of the spins remain in a frozen state for 
small system sizes. Assuming then that the nematic state 
in a finite system is characterized by this residual corre- 
lation, we can define a normalized correlation function 
C'{t) as 



C'{t) 



{C{t) — Cplateau) 
(C'(l) ~ Cplateau) 



(5) 



where C'{t) — characterizes the equihbrium nematic 
state (either stable or metastable) . In Fig[T6] we see 
the typical behavior of C'(t) for different temperatures 
from above T2 down to temperatures well below Ti . Also 
C'{t) exhibit different behaviors above and below T2. For 
temperatures above T2 the curves are very well fitted (see 
FigfT6| by a modified stretched exponential function of 
the type 



f{t) = A re- 



-(t/rr 



(6) 



This fact is peculiar of the tetragonal Hquid state: it 
does not exhibit exponential decay, even at relatively 
high temperatures. This behavior is observed also for 
S ~ 1. For temperatures below Ti, C'{t) exhibit two 
different regimes, as can be appreciated more clearly in 
FigflTl At short times C'{t) is again well fitted by a 
function of the type Eq. ^ . For times longer that some 
crossover time, C'{t) enters into a pure stretched expo- 
nential regime, i.e., the best fit is obtained by a function 
of the type 



/(t)=Ae-(*/^)^ 
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FIG. 16: (Color online) Normalized correlation function 
Eq.im as a function of t for 6 = 2 adn L = 36, after a quasi- 
static cooling at different temperatures. The continuous lines 
for temperatures above T2 correspond to fittings using Eq.®. 
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FIG. 17: Normalized correlation function Eq.([5]) as a func- 
tion of t for 5 = 2, after a quasi-static cooling at T = 0.78. 
The continuous and dashed lines correspond to fittings using 
Eqs.© and ([7} respectively. 



Finally, we analyzed the temperature dependence of 
the relaxation time t. Since there is too much uncer- 
tainty in fitting a curve with three parameters, instead 
of considering r we defined another correlation time r* 
as the time at which C"(t*) = 0.1. In Fig [18] we see 
an Arrhenius plot of r* for different system sizes. The 
dashed lines in the figure correspond to fits using a Vogel- 
Fulcher-Tamman (VFT) form 



r = To cxp 



A 



T-Tn 



(8) 



(7) 



which is usually associated with the behavior of relax- 
ation times in a fragile glass-former, according to Angell's 
classification^. The fitting values of Tq for the different 
system sizes are shown in Table [ill We also show the 



n 



L 


To 


K 


36 


0.41 ±0.025 


1.20 ±0.4 


48 


0.40 ±0.015 


1.25 ±0.3 


60 


0.44 ± 0.075 


2.10 ±0.8 



TABLE 11: Fitting values of To and the fragility parameter 
K = To/ A for the relaxation time r* as a function of the 
supercooling temperature using a VFT form. 
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FIG. 18: (Color online) Relaxation time r* vs. l/T for S = 2 
and different system sizes. The dashed lines correspond to 
fittings using VFT form. The vertical line corresponds to 
1/Ti 



fragility parameter K = Tq / A, which show very high val- 
ues around « 1. It is worth mentioning that the data 
of FigdHl can also be fitted with a power law of the type 
T* = A\T — Tq\~^. However, the best fitting is obtained 
with the VFT form Eq.®. 

Notice the strong dependence of r* with L for the range 
of sizes here considered. In Fig[T9]we show t* as a func- 
tion of L for a fixed temperature in the metastable region. 
We see that r* becomes independent of L for sizes larger 
than i ^ 80, which coincides with the crossover value 
observed for the plateau, suggesting the same origin for 
both finite size effects. 



IV. DISCUSSION AND CONCLUSIONS 

The main result of this work is the observation of inter- 
play between nucleation and coarsening which appears to 
be a novel feature. In other words, coarsening happens 
usually as a competition between domains of different 
stable phases. In the present case both direct coarsen- 
ing of the stripe phase and nucleation of the same phase 
on a nematic background can happen with a probabil- 
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FIG. 19: Relaxation time t* vs. L for 5 = 2 and T = 0.7. 



ity that depends on the final temperature. Nucleation 
processes are associated with a strong metastability of 
the nematic phase. If the final temperature T/ < Ti 
is not too far below the stripe-nematic transition tem- 
perature the probability of nucleation is high: during 
the quench the system initially relaxes to the nematic 
phase, which is metastable below Ti until finally it de- 
cays through nucleation to the equilibrium stripe phase. 
For final temperatures deep in the stripe phase simple 
coarsening processes between stripe states with different 
orientation have larger probability of occurrence than nu- 
cleation. Although no detailed theory is available for 
predicting this behavior in the whole temperature range, 
classical arguments assuming homogeneous nucleation, 
coarsening with a non-conserved order parameter and fi- 
nite size scaling give a reasonably good interpretation 
of all the observed phenomenology. The origin of this 
interplay may rely in the fact that the stable and the 
metastable phases share the same orientational symme- 
try, so we can have domains with stripes (or pieces of 
stripes) that can be oriented along the same directions of 
the square lattice substrate in both phases. 

Quasi-equilibrium autocorrelations in the stable and 
metastable nematic regions show a plateau whose nature 
is not completely clear. For temperatures where the ne- 
matic phase is only metastable, this plateau appears in a 
finite time window until a final decay for times t > Tnuci ■ 
For temperatures in the equilibrium nematic phase the 
behavior is more subtle. Although in a first approxima- 
tion one could consider it to be a distinct characteristic 
of the nematic phase, in which a fraction of the spins re- 
main practically frozen, we got indications that for the 
larger sizes the height of the plateau begins to diminish, 
signaUing a possible finite size origin. One possibility is 
that the fraction of frozen spins in the nematic phase di- 
minish with system size, in which case the plateau will 
eventually disappear for sizes large enough. More work is 
needed to clarify the behavior of equihbrium relaxations 
in large samples. 

The present results also confirm the existence of the 
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intermediate nematic phase, as observed in equilibrium 
simulations iiii^ii^; all the dynamical observations are in 
agreement with the reported equilibrium results. In fact, 
a nematic phase with quasi-long-range order is expected 
to be generic in isotropic systems with competing inter- 
actions^. Recent analytic results on a continous version 
of the present system show the presence of an isotropic- 
nematic transition which is in the Kosterlitz-Thouless 
universality class, while the stripe phase with true po- 
sitional long range order turns out to be unstable for 
continous and isotropic systems in d = 32i in the thermo- 
dynamic limit, at variance with the discrete model stud- 
ied here. Nevertheless, it is important to note that, al- 
though truly unstable in the thermodynamic limit, stripe 
domains grow already in the paramagnetic phase with 
decreasing temperature and may be relevant for ob- 
servations in finte size systems, like in simulations for 
example. In the present case, the stability of the low 
temperature stripe phase may be due to the symmetry 
breaking effect of the lattice and the sharp nature of the 
Ising domain walls. A unified interpretation of the nature 
of the transitions in the discrete and continous models is 
still lacking. 

A final important point concerns the strong increase 
observed in the stationary relaxation time r* in the 



metastable nematic phase. Although a spinodal insta- 
bility cannot be excluded, the stretched exponential re- 
laxation and the super Arrhenius behavior of r* (T) point 
to an increasing cooperative and glassy behaviour of the 
dynamics. Indeed, a very similar behavior in a related 
model, namely the Coulomb frustrated model in 3d, has 
been interpreted as an evidence of fragile glass forming 
phenomenologj^. Recent experimental results on Fe on 
Cu(lOO) ultrathin films^ have suggested the possibility of 
a stripe-liquid to stripe-glass transition in the vicinity of 
the stripe ordering temperature. Stripe-glassy behaviour 
can be interpreted within the Frustration Limited Do- 
main Theory21 or the Random First Order Transition 
— scenarios. We have not attempted a quantitative 
comparison with either theory, which remains as subject 
of future work. 
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